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Context. Thanks to their high quality, new and upcoming asteroseismic observations - with CoRoT, Kepler, and from the ground... - 
can benefit from the experience gained with helioseismology. 

Aims. We focus in this paper on solar-like oscillations, for which the inclination of the rotation axis is unknown. We present a 
theoretical study of the errors of p-mode parameters determined by means of a maximum-likelihood estimator, and we also analyze 
correlations and biases. 

Methods. We have used different, complementary approaches: we have performed either semi-analytical computation of the Hessian 
matrix, fitting of single mean profiles, or Monte Carlo simulations. 

Results. We give first analytical approximations for the errors of frequency, inclination and rotational splitting. The determination of 
the inclination is very challenging for the common case of slow rotators (like the Sun), making difficult the determination of a reliable 
rotational splitting. Moreover, due to the numerous correlations, biases - more or less significant - can appear in the determination of 
various parameters in the case of bad inclination fittings, especially when a locking at 90° occurs. This issue concerning inclination 
locking is also discussed. Nevertheless, the central frequency and some derived parameters such as the total power of the mode are 
free of such biases. 

Key words. Stars: oscillations - Methods: data analysis - Methods: statistical 



1. Introduction 

A new era of the asteroseismic observations is now opening with 
CoRoT (Baglin et al. 2006), and will be pursued in the near fu- 
ture with Kepler (Borucki et al. 2007). As these missions will 
provide long and uninterrupted time series of intensity measure- 
ments, the quality of these data will become closer to those of 
helioseismology, although with a higher level of noise expected. 
In such a context, asteroseismology can fruitfully inherit the 
methods and techniques developed for helioseismology, espe- 
cially for so-called 'Sun-as-a-star' observations, such as those 
provided by the space instruments GOLF (Gabriel et al. 1995) 
or VIRGO (Frohlich et al. 1995) onboard the SoHO spacecraft 
or by ground-based networks like BiSON (Chaplin et al. 1996). 

We discuss in this paper the usual way to extract information 
on stochastically excited acoustic (p) modes in helioseismology. 
P modes can be described by several parameters (frequency, life- 
time, amplitude... ) which can be derived by fitting the stellar 
oscillation power spectrum with a maximum-likelihood estima- 
tion. As for any measurement, we need to associate with each 
derived parameter a correct error bar, that is essential to estimate 
the significance of the measurement and thus to be able to do a 
reasonable interpretation. We present here theoretical results on 
the derivation of errors of p-mode parameters and their correla- 
tions. This work generalizes some results already known in he- 
lioseismology (Libbrecht 1992; Toutain & Appourchaux 1994) 
by adding the inclination of the stellar rotation axis, i, as an extra 
free parameter. 



Fitting methods with a free inclination angle have first been 
analyzed by Gizon & Solanki (2003), then by Ballot et al. 
(2006). In the present paper, we develop a more complete version 
of a preliminary work (Ballot et al. 2008) on analytical formu- 
lations of error bars for several parameters, and we have com- 
pleted it with an analysis of correlations and biases. We have 
especially considered here the very common cases of slow ro- 
tators, for which the mode linewidth is greater than (or similar 
to) the rotational splitting, giving rise to a blending of multiplet 
components. 

In § 2, we define the model assumed for p modes and the 
fitting method. Section 3 deals with error bars obtained for the 
mode frequencies, the splittings and the inclination. The corre- 
lation between different parameters is studied in § 4 and some 
biases of the method are analyzed in § 5. Finally we briefly dis- 
cuss the locking of the inclination determination at 90° which 
often appears during fitting (§ 6). 

2. Models and Methods 

2. 1 . Modeling the power spectrum 

Stellar acoustic eigenmodes are characterized by their degree I, 
their azimuthal order m and their radial order n. In this study, 
we treat p modes according to the solar paradigm. Modes are 
modelled as stochastically excited and intrinsically damped har- 
monic oscillators (see Kumar et al. 1988). In that case, the power 
spectrum - obtained by computing the discrete Fourier transform 
of an evenly sampled timeseries - of such modes is distributed 
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around a mean Lorentzian profile with an exponential probabil- 
ity distribution. A Lorentzian profile is defined as 



L(A, vo, T; v) 



1 + 



(i) 



where A is the mode height, T = (tit) -1 the mode linewidth, 
linked to the damping time r, and vo the mode frequency. For our 
examples we have considered the value F = 1 fiHz in the whole 
study, that corresponds to a lifetime around 3^4 days. This is 
a typical observed value of solar p-mode linewidths in a broad 
range around 2500-3000 fiHz (e.g. Garcia et al. 2004). 

Since stellar rotation lifts the azimuthal degeneracy of eigen- 
modes, the power spectrum of a mode with a degree I is a multi- 
plet with 21+1 components. Extrapolating from the solar case, 
we assume equipartition of energy between the components in a 
multiplet, and we define a multiplet as a symmetric profile: 



Mi(A, v , T, v s , i; v) = / \ a Lm (i)L(A, v + mv s , F; v), 



(2) 



m--l 



where v s is the rotational splitting and the inclination ; is the 
angle between the rotation axis and the line of sight. 

We have used here the approximation of Ledoux (1951) for 
an uniform rotation: the frequency shift of the m-component rel- 
ative to the central one is v,„ - vq = mv s \ moreover, we consider 
v s « Q./(2n), that is the asymptotic regime at high order n (Q is 
the stellar rotation rate). This approximation is valid when rota- 
tion is sufficiently slow and there is neither a strong differential 
rotation nor a strong magnetic activity. 

Next, the inclination acts only on «/,„,(/), the amplitude ratios 
of components inside a multiplet, which satisfies the relationship 
Tj m a i.m = 1- Thus the total power of a multiplet is always P = 
% AT. The amplitude ratios are written 



|/-/n|! 

ai,m(i) = 7, rr(f/,/«(cos 0) 

\l + m\\ 



(3) 



where Pi m are the associated Legendre functions (see Gizon & 
Solanki 2003). To derive this expression, when the flux is in- 
tegrated over the full stellar disc, we need to assume that the 
weighting function, which gives the contribution of a point on 
the disc to the integral, is a function of the distance to the disc 
center only. This is correct for intensity measurements, because 
the weighting function is mainly linked to the limb-darkening. 
However for velocity measurements, we can observe departures 
from this law, since the rotation of the star induces an asymmetry 
in the velocity field. 

Last, we have also assumed that a mode is not correlated with 
any other modes or with the convective background noise. Doing 
so we neglect in the present study any possible asymmetry of 
the Lorentzian profiles (Nigam et al. 1998). According to what 
we have learned from helioseismology, neglecting asymmetries 
could introduce systematic errors (i.e. biases) in mode frequency 
determination. These errors are of the order of 0. 1 juHz in the so- 
lar case (Toutain et al. 1997). This is on par with the statistical 
error (i.e. the standard deviation) of the frequency for time series 
of several months. For longer time series (few years), the asym- 
metry should be included in the fitted profiles to avoid systematic 
errors in mode frequencies. 

To summarize, we have considered that the star is observed 
in intensity, that the mode excitation mechanisms are close to 
those of the Sun, and that the star does not rotate too rapidly (a 
few times the solar rate). 



2.2. Maximum-likelihood estimator 

P-mode parameters are fitted with a classical maximum- 
likelihood estimation technique and the associated error and cor- 
relation are estimated by inverting the Hessian matrix (Toutain 
& Appourchaux 1994; Appourchaux et al. 1998). In practice, in- 
stead of maximizing the likelihood, we minimize the negative 
logarithm of the likelihood function, which yields for a random 
exponential noise: 



lnXW = Y j \nS(A;v k ) + 



S{A-v k ) 



(4) 



where S is the model of the spectrum measured {S k}k=\.N at fre- 
quencies {v/J, and A = (Ai, . . . , A p ) is the set of p parameters to 
be adjusted. We denote hereafter Aj an estimation of Aj. 

As the observed intensity fluctuations are integrated over the 
whole stellar disc, only low-degree modes are visible. A quick 
estimation of their visibility indicates that we will mainly be able 
to detect only modes with I < 2 and perhaps a few / = 3 modes. 

We have therefore considered the following fitting cases: 1) 
we fit a mode / = 1 alone; S is described by 6 parameters: A, vq, 
T, i, v s , and a additive background B assumed to be flat within 
the fitting window. In practice, we fit the logarithm of some of 
the parameters: a - In A, y — lnT and b = \r\ B. 2) We fit a pair 
of modes / = and 2; assuming a common linewidth for both 
modes S is described by 8 parameters; 3) We fit a sequence of 
modes I = 0, 2 and 1 ; S is described by 1 1 parameters, assuming 
the splitting is the same for the consecutive / = 1 and 2 modes. 

The standard deviations (error bars) crj associated with Aj 
are estimated with the covariance matrix C, computed by invert- 
ing the Hessian matrix H (C = H~'). The diagonal elements 
of C give the errors cjj = a 1 -, while the non-diagonal elements 
give the covariances cy = cry = pij<Ti<Tj (py are the correlation 
coefficients). 

The terms of the Hessian are: 



hij = 



d 2 i(A) 



dAidAj 



(5) 



Following Libbrecht (1992) and Toutain & Appourchaux 
(1994), we define a theoretical Hessian corresponding to the av- 
erage of a large number of realizations: 



1 



dS dS 



k=\ 



S 2 (A;v k ) dAidAj' 



(6) 



If the frequency bin is much smaller than the mode linewidth, 
we can approximate the Hessian by the integral: 



1 



dS dS 



S 2 (A;v) dA^Aj 



dv, 



(7) 



where T is the observation duration. 

In the next sections, results concerning errors, correlations 
and biases are obtained with semi-analytical computations of 
Eqs. 6 and 7, by fitting the mode profile model as in Toutain 
et al. (2005), or with Monte Carlo simulations. 

3. Theoretical error bars 

3. 1. Error of the central frequency vq 

We generalize in this section the results of Libbrecht (1992) and 
Toutain & Appourchaux (1994), inferred for / = and for mul- 
tiplets in the solar configuration (; = 90°). 
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Fig. 1. Frequency error computed with Eq. (9) for modes / = 1 (left) and 1 = 2 (right) as a function of the angle i, for different reduced splittings 
x s = 0.4 (dots), 0.8 (solid line), 1.6 (dashes), 2.4 (dot-dash), 4.0 (dot-dot-dot-dash), 8.0 (long dashes); that corresponds, for T = 1 yuHz, to Q. =0.5, 
1, 2, 3, 5, 10 Cl . Thick long dashes indicate the limit when x s » 1, computed with Eq. (13). Here /? = 1/20 and T = 6 months. 



Considering an isolated multiplet S — B + Mi, we can easily 
verify that h Vo j(= /^> ) = except for Aj = vo. Since S is even in 
the variable v-Vo, the derivatives d Vo S = —d v S are odd. However 
we can verify that is even relative to the variable v - vo for 
all of the other parameters Aj — B,A,F,v s and i. This is easy 
to demonstrate by noting that any modification of one or several 
of these parameters preserves the symmetry of the profile with 
respect to v-Vo. With the same argument, we easily demonstrate 
that, after any variable change Aj = f(B,A, F, v s , i), d^S is still 
even. Thus, in Eq. 7, the function under the integral is odd, and 
the integral vanishes. We deduce that 



and 



cr Vaj = VAj + v . 



(8) 



The lack of correlation of vo with all other parameter is well es- 
tablished in Monte Carlo (MC) simulations. When a pair (for 
instance I = and 2) is fitted, this is still correct while the 



modes are well separated: given A^ and X£' two parameters of 

the considered I = and I = 2 modes, then dmSd^S ~ Vv, 

j * 

since either 3,<o>S or d A a)S vanishes, as modes are sufficiently 

j * 

far apart. However, when modes are blended, crosstalk can oc- 
cur and the central frequency can become correlated with other 
parameters, especially with the central frequency of the neighbor 
mode. Moreover, if the mode profiles are not symmetric (asym- 
metries in the Lorentzian or in the splittings), the mode frequen- 
cies can also become correlated with the others parameters, but 
a priori we expect negligibly small effects. 
From Eq. 7, we derive the error for vo: 

or 2 VB = ~mx s ,i), (9) 

with B — B/A the noise-to-signal ratio and x s = 2v s /r the re- 
duced splitting. For I — 0, we find the formula of Libbrecht 
(1992): 



l(2) 



/oOSonly) = #TT(#TT+ jfi) . 

For / > 0, the simplest form of // is its integral form: 

2m ai,m{i)(x + mx s )L 2 r (x + mx s ) 



(10) 



mx s ,i) = - 

where 

L r (x) 







P + 2m ai Jn (i)L r (x + mx s ) 



dx 



is the reduced Lorentzian. 

Figure 1 shows the evolution of cr Vo with ; and x s for I - I and 
2 modes. As expected, /;(/?, x s , i) approches fo(B), as ; or x s ap- 
proches zero. We can also notice that, when x s <k 1, the splitting 
acts as an extra width, increasing the error for vo. Furthermore 
we clearly see that, depending of the values of i and x s , the error 
bars for vo can vary within a factor of 2. When x s » 1, i.e. for 
large rotation rates, the components are well-separated, the fit- 
ting gives exactly the same result as when each ;7?-component is 
independently fitted and the central frequency is computed with 
a weighted average: vo = 2m 'u"', 2 / 2m °~v 2 - The associated 
error is then 



^ = — L 

m=-l m=—l 



fo 



ai, m (i) 



(13) 



This value, plotted in Fig. 1 as a thick dashed line, gives the 
lowest limit for the error at a given angle, linewidth and S/N 
ratio. 



3.2. Error of the angle i and the splitting v s 

We focus now on the error of the inclination and the splitting, 
or other related variables. Ballot et al. (2006) have used the pair 
of parameters (i, v*) instead of (i, v s ), where v* = v s sin i is the 
projected splitting. The authors did not see any difference in the 
determination of i using one set of variables or the other. We 
analytically derive this result here, and we study more generally 
what happens when we fit other sets of parameters, for instance 
(sin 2 i, v^) as further mentioned in § 6. 

Let us denote by (p, q) a new pair of parameters such i = 
fi(p) and v s = fs(p, <?)■ We consider a simplified Hessian matrix 
2x2 with elements (see Eq. 6): 



(1 



1 dS dS 
The errors for i and v s are 



1 + x 



(12) 



-2 



h 2 

1 _ (£ 

h. ' 



(14) 

(15) 
(16) 

(17) 
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Fig. 2. Error for the angle i, the splitting v„ and the projected splitting v*, as a function of ;', for modes 1=1 (left) and 1 = 2 (right), for different 
values of x„ (see caption of Fig. 1). Black lines show errors obtained with the simplified Hessian, thin red lines are obtained by inverting the full 
Hessian. Thick long dashes in the middle panels indicate the limit of the error on v, when x s » 1, computed with Eq. (22). Here fi = 1/20 and 
T = 6 months. 



-2 



hi 
hi 



(18) 



We denote by G the simplified Hessian matrix for the new pa- 
rameters. Its elements are g p , g q and g pq . Using rules of (partial) 
derivatives, we show that: 



hih s - h is 



dp 
di 



and 



<df s dfi 




dq dp ] y hi 
Specifically for p = i and q = v*, we find: 



ov = o~ v sin i -i / 1 H , 

V tan z i hi tan ; h t 



-2 



v s h is 



(19) 



(20) 



(21) 



First we conclude that the error for ; does not depend on 
the choice made for the other variables: v s , v*, or another com- 
bination. MC simulations and all of our other tests confirm this. 



Second, using cos ; or sin i does not modify the error: we retrieve 
natural relationships such as cr sin i i = | sin2/|<x,. 

Figure 2 shows the errors we derive for the angle i, the split- 
ting v s and the projected splitting v* for 6-month long obser- 
vations. We compare the results from a simplified Hessian to 
those obtained by fitting directly the profile model (see method 
in Toutain et al. 2005). For the angle, both computations give 
very similar results. Estimated errors clearly demonstrate the dif- 
ficulty in deriving a reliable determination of the angle ; from 
one single mode when the rotation is less than 2Qq: the uncer- 
tainty covers almost the entire possible range. / = 2 modes pro- 
vide slightly better estimates of ; than / = 1 modes because they 
have more components and the displacement of the / = 2 modes 
is twice that of the / = 1 modes. These error estimations are in 
agreement with those of Gizon & Solanki (2003), who have also 
numerically derived errors from the Hessian. Dividing by V8, 
the values of the errors in Figure 2 would give the correspond- 
ing errors for 4 years of observation (scales as the square root 
of the ratio of the observing times). Unless several modes are 
used to determine the inclination and the rotation, 4 years is still 
not sufficient to get enough accuracy when Q % Oq. However, 
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if we measure 10 independent modes (for instance, 5 I — 1 and 
5 1 = 2), by averaging the results, we can expect to reduce the 
error to ~ 15° for Q = Q0. 

Concerning the error for the splitting v s , our simplified anal- 
ysis also gives a result in good agreement with the complete nu- 
merical computation. We notice as well the dramatic increase of 
the errors when v s decreases. As previously done for the central 
frequency, we can also derive here the limit of the splitting er- 
ror when x s » 1 and that all of the components can be fitted 
independently. The resulting error for v s in this case is given by: 



cr 



?-Z' 

m——l 



8ttT 



m cr„ = 



/o 



fl;. m (0 



(22) 



Interestingly enough it gives a lower bound for cr Vj for given i, 
T, T and /3. 

The lower plots in Fig. 2 show the error of the projected split- 
ting v*. One sees that our simplified calculation is slightly more 
crude in this case. In this simplified approach, the correlations 
of i, v s and v* with the other parameters are neglected. They are 
actually non-negligible, especially for the linewidth F. For slow 
rotation, v* is a bit more correlated with the linewidth F than 
v s (see § 4), that mainly explains why the discrepancy between 
simplified and full Hessian computations is larger for v* than for 
v s . Nevertheless, the most remarkable fact is the smaller scatter 
of the curves compared to the previous plots: as suggested by 
Ballot et al. (2006), this computation indeed shows that cr y > is 
noticeably less sensitive than cr Vj to the value of v s . 

4. Correlations 

We discuss in this section the correlations between the different 
parameters of a mode and the effects of the correlations on data 
analysis. 

Let us assume as an example for the remaining of this paper 
a star characterized by i = 45°, v s = 0.8 pHz {2Q.q), F = 1 //Hz, 
and /3 = 1/20 for / = 1 modes and 1/10 for I = 2. This cor- 
responds to signal-to-noise ratios about five times smaller than 
those observed for the Sun with VIRGO. We consider here a 
single / = 1 multiplet observed for 4 years - that corresponds 
to Kepler long observation runs - and perform a MC simulation 
with 10000 realizations. The theoretical errors and correlations 
have been deduced from the theoretical Hessian. Figure 3 shows 
both MC and theoretical results and allow us to verify the con- 
sistency of both computations for several pairs of parameters. 
MC results are shown as clouds of points. From the theoretical 
Hessian, we have derived the error bars, but also the ellipses of 
errors and the regression lines to make visible the correlation. 
We recall that an ellipse of errors is, for a 2D normal distribu- 
tion, an isoprobability curve given by the equation Q(A,, Aj) = k 2 
where Q is the quadratic 



Q(Ai,Aj) = 
1 



1 



■Pij 



(Ai - m) 1 2p u (Ai - m)(Aj - m) (Aj - hj) 2 



cr 2 



trier 



(23) 



where yu,, cr, are the mean values and the standard deviations and 
Pij the correlation coefficient (see §2.2). The number k is a real 
constant - generally integer - and defines the confidence level: 
the probability to get a point inside the ellipse is 1 - exp(-^A: 2 ). 

As for the Sun, the pair of parameters (A, T) are strongly cor- 
related such AY ss const (not plotted). We have found a correla- 
tion coefficient p ay = -0.91 for our example. This coefficient is 



almost insensitive to the values of the inclination and the split- 
ting. By exploring different values of i and v s , we always find 
a correlation with an absolute value larger than 0.8. This is, of 
course, independent of choosing either v s or v* as a free param- 
eter. 

The second strong correlation we observe is between i and 
v s (Fig. 3a). We recover the result of Ballot et al. (2006): for low 
rotation rate, (i, v s ) are correlated such that v s sin i » const. In 
our example, the correlation coefficient p; Vj = -0.68, and |p lv J 
is even greater than 0.9 for Q = Q . By using v* instead of v s 
(Fig. 3b), the correlation decreases: p,y = 0.21 for our example, 
and we have obtained p iv j ~ -0.1 for Q. = Q Q . It is important 
to recall that it is true at low rotation rate only; when x s » 1, 
the situation is opposite, ; is more correlated with v* s than with 
Vj. Figure 4 illustrates, for a I = 1 mode, the change of p IVj from 
-1 to as v s increases. We get similar results for / = 2. The 
correlation coefficients are deduced from both the full theoretical 
Hessian and the simplified one proposed in §3.2. We note the 
limitation of the latter calculation for intermediate values of v s : 
for these configuration, that gives the good order of magnitude 
for pi Vs , but an incorrect dependency on i. 



0.4 7 

0.2 : 

0.0 - 
-0.2 - 
-0.4 - 
-0.6 
-0.8 



1=1 



-l.Ot 







20 



40 



60 



80 



i(°) 



Fig. 4. Correlation coefficient between i and v, for an / = 1 mode as 
a function of (, for different values of x s (see caption of Fig. 1). Black 
lines show errors obtained with the simplified Hessian, thin red lines are 
obtained by inverting the full Hessian. 

Figures 3c and 3d illustrate the correlation of T with v s and 
v*. We can also see it by looking at the color gradient, which 
appears on the upper panel. This correlation is significant and 
is larger with v* than with v s . That is understandable: the fitting 
technique has problems distinguishing between the broadening 
due to splitting and that due to the natural linewidth. We have 
verified that the correlation decreases strongly when x s » 1 . 

Let us finally notice the dense group of points with i * 90° 
obtained with the MC simulation (Figs. 3a and b). The phe- 
nomenon is discussed in § 6. These peculiar fits are indicated 
by red dots on Figs. 3c and 3d and we clearly note that, due to 
correlations, this set of points is shifted - this indicates a notice- 
able bias especially on v s . We notice also these points are orga- 
nized along the regression line at constant v s , which indicates 
that, when i is locked at 90°, v s is almost blocked around an 
underestimated values; that is another consequence of the (i,v s ) 
correlation. 

Ballot et al. (2006) suggest using v*, instead of v s , at low ro- 
tation rate, to improve some averages and avoid some effects of 
correlations: for instance to estimate amean splitting, (v*)/(sin/) 
can be a better estimator than (v s ). We have verified this point 
with a MC simulation of our example. We have performed 2000 
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Fig. 3. Top: correlations between the angle i and the splitting v, (a) or the projected splitting v* (b). Dashed line indicates v* = const (a) or 
v, = const (b). Point color indicates the fitted value of linewidth y. The color scale is saturated at ±cr y (<r y = 0.108). Bottom: correlations between 
the linewidth y and the splitting v, (c) or the projected splitting v* (d). Red dots indicate fits where i > 89°. For each plot, lcr error bars, ellipses of 
errors (with k = 1, 2 and 3; see text) and regression lines are deduced from the mean Hessian, assuming normal distributions of results. The value 
of the correlation coefficient p deduced from the mean Hessian is also shown. 



realizations of 10 modes (5 / = 1 modes and 5 / = 2 modes). 
The distribution of (v*)/(sin/) compared (v s ) exhibits a smaller 
dispersion (0.10 against 0.15 pMz in our case) with a reduced tail 
at high v s . 

One must also take into account the existing correlations 
when one computes derived variables such as the total power 
of a mode P. We denote p — In P. A naive estimation of cr 2 is 
a 2 + cr 2 » 0.47 2 (for our example with T — 6 months). However 
due to the correlation, we have to consider the covariance to re- 
cover a correct error: cr 2 — cr 2 + cr 2 + 2cr a7 ~ 0.15 2 (this value 
is in agreement with those found with a MC simulation, see §5 
Fig. 5). This point seems obvious but is frequently forgotten. 

Lastly, when it is possible to determine the angle i with other 
techniques, one can fix it and thus one stands in a position simi- 
lar to the helioseismology, which is more confortable. However, 
we have to keep in mind that, due to the complex correlations 
between all of the parameters, an angle i assigned to a wrong 
value introduces biases in the determination of other parameters 
and so introduces systematics which should be estimated on a 
case by case basis. 

5. Biases for A and cr 

We know that asymptotically, for long observation times, our 
estimators of parameters Aj and their associated errors &j are 
non-biased (e.g. Appourchaux et al. 1998). However, because 
observation durations are finite, biases can appear both on Aj 
and &j. To study this point, we have used the results of the MC 
simulation described in the previous section. Asymptotically, for 



T — > oo, due to the Central Limit Theorem (CLT), fitted parame- 
ters follow normal laws; thus (Aj-Aj)/d m j asymptotically follows 
a standard normal law. Figure 5 shows, for our simulated case, 
the distribution of the ratio between real errors Aj - Aj and esti- 
mated error bars cr ; for different parameters and for two different 
observation durations: 6 months (close to CoRoT long runs) and 
4 years (Kepler long runs). These plots allow us to see potential 
biases on parameter determinations or over/underestimations of 
errors, by comparing the position and the shape of the distribu- 
tion relative to the standard normal law. 

First, the logarithmic linewidth y is slightly biased after 6 
months of observations: fitting tends to underestimate its value. 
However, the distribution spread is close to the expected one, 
that indicates correct estimations of errors <x r . The distribution 
of logarithmic heights a (not plotted) is very similar to those of y, 
and in particular we find a similar small bias for short times, but 
in the opposite sense. As a consequence of these opposite biases, 
the sum of both, i.e. the mode power p, is almost non-biased, 
even for 6-month observations and its distribution is symmetric 
and almost normal. If we turn now to the central frequency, vo, 
we notice no bias for the fitted values vo, but for 6-month obser- 
vations, the errors are generally underestimated, as shown by the 
extended wings of the distribution. 

Six months are generally too short to verify the CLT, espe- 
cially for splitting - projected or not - and angle. The distribu- 
tion of v s is non-gaussian and is moreover highly spread, that 
indicates cr Vt are significantly underestimated. The distributions 
of i and v* are less spread than the previous one, but exhibit clear 
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Fig. 5. Distribution of real errors, Aj - Aj, normalized by the estimated 
errors <x, for different parameters (from top to bottom: linewidth, cen- 
tral frequency, splitting, projected splitting, angle and mode power), for 
two different durations: T = 6 months (left) and 4 years (right). On 
each plot, the filled histogram is the distribution when only realizations 
with i > 89° are taken into account, dashed line is the standard nor- 
mal distribution, and the mean value (cry) is the error deduced from the 
theoretical Hessian, quoted in fjHz for v , v, and v*, and in degrees for 



asymmetries, due to biases for both estimated parameters and er- 
rors. 



Table 1. Fraction of fitting locked at i = 90° for different observation 
times and different sequences of fitted modes 
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After 4 years, all of the distributions are close to normal. 
Nevertheless, the angle distribution is still asymmetrical and ex- 
hibits a slight bias. The distributions of splitting and projected 
splitting also show a very small asymmetry but become very 
close to the standard gaussian, that is a noticeable change mostly 
for v s . 

Let us focus now on the realizations for which fits have been 
locked at 90°. First, we note that the number of such cases is 
reduced for T = 4 years compared to 6 months. Due to correla- 
tions, when the fitting converges to ; = 90°, all of the parameters 
are expected to be biased; there are few exceptions: the central 
frequency, because it is not correlated with the other parameters 
(§3.1), and the mode power, a derived parameter which appears 
to be almost non-correlated with the angle. Another manner to 
understand this fact is that, since i only acts on the shape of the 
multiplet, it does not modify its integral (i.e. the power P) or 
its center of gravity (i.e. vo). All other parameters are affected, 
and the bias becomes significantly large for v, (greater than 5<x) 
when T — 4 years, which is visible as a bump on the left of 
the plot. It has already been mentioned in §4 and is visible on 
Fig. 3c as a set of crosses shifted to the left toward low split- 
ting by about 3<x. These two results seem in contradiction at first 
sight, but <x has a different meaning in Fig. 3c and Fig. 5. From 
Fig. 3c we learn that, when locking occurs, the mean bias for 
v s is 3 times the mean standard deviation since red dots cluster 
around the outer ellipse of errors corresponding to the 3<x level, 
whereas Fig. 5 shows that, in such situations, the mean bias for 
v s is 5 times the error bar estimated for the given case. This in- 
dicates that, when locking occurs, the error bars for v s also are 
underestimated. 

6. Fit locking at i = 90° 

Previously, we have seen that a non-negligible proportion of fit- 
tings converge toward a solution with i = 90°. This phenomenon 
has been first observed by Gizon & Solanki (2003). We summa- 
rize in Table 1 the fraction of such fits for different situations. 
Once again, we have considered here our sample star, for several 
observation durations from 3 months to 4 years. We have fitted 
single / = 1 modes, pairs of modes / = and 2, and also se- 
quences of three modes / = 0, 2 and 1 . I — 1 modes have been 
fitted using (i, v s ) or (i, v*) as free parameters, and, in one case, 
all of the other parameters have been fixed to their exact values. 

We can draw several conclusions. As we expect, the frac- 
tion decreases when T increases; thus, for Kepler-like observa- 
tions (4 years), there is almost no problem. Second, the results 
are almost the same by fitting either (/, v s ) or (i, v*), there is no 
improvement by using v*. We notice also that even when we 
impose the input A, F and v , the part of bad fittings is not sig- 
nificatively reduced. This indicates the problem is intrinsic and 
not just a convergence failure as explained below. Next, since 
they have more components, / = 2 modes give slightly better re- 
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suits, though their S/N ratio is lower. However, for I = 2, another 
problem appears: in numerous configurations, the fitting hardly 
distinguishes between a low v s with a high i, or a lower i with 
a v s two times larger, due to a confusion between m = +1 and 
+2 components (See also Gizon & Solanki 2003). Fitting simul- 
taneously / = 1 and 2 modes allows us to solve this ambiguity. 
However, although the best results are obtained with sequences 
of three modes, there are still almost one bad fitting out of ten, 
for 6-month observations (CoRoT-like observations) - that is not 
negligible. 
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Fig. 6. Maps of likelihood t = -lnX in the plane (/, v.,) for a / = 1 
mode. Top: map of the mean likelihood. Bottom: an example of realiza- 
tion where the locking phenomeoccurs. Black corresponds to the min- 
imum of t, white to higher values. Fine white or black dotted lines are 
contours of constant likelihood. Crosses indicate the input values of i 
and v, in the spectum; dashed line indicates y* = const. Solid lines in- 
dicate the minima of C along i at fixed v„ i.e. points where d\i = 
and dj£ > 0; dot-dashed lines indicate the minima of I at fixed i; thick 
lines mean that this is a global minimum, fine lines a local one. On the 
bottom plot, the cross with errors bars indicates the result of the fit. 



To understand the origin of the fit locking, we have plot- 
ted on Fig. 6 {top) the mean likelihood in the plane (i, v s ) for a 
/ = 1 mode. The plot is symmetric with respect to 7r/2 (and n- 
periodic). As expected, minima correspond to the values of (i, v„) 
used to model the spectrum. Since the likelihood, I, is even in 
i - n /2, a valley along ; = n/2 appears at low splitting, before 
a fork appears at higher splittings. Obviously the (black) region 
around the minimum extends up to the limit of 90°, which means 
that it is quite likely that for a noise realization the minimum of 
the likelihood will end up close to that limit. It may even end 
up beyond that limit and in this case because of the periodicity 
of the inclination there is obviously no inclination which could 
allow the fit to reach such a minimum. Figure 6 bottom shows an 
example of such a situation with a global minimum at 90°(and 
a secondary minimum around 20°, with interestingly the same 



value of v* since it is also on the dash line). The fit has then 
been blocked by the hard limit at 90° and has stopped there pro- 
ducing this effect of locking at 90°. A similar effect is described 
in Chaplin et al. (2008) for a problem of fitting solar-cycle fre- 
quency shift or in an ongoing work by the same authors about 
the zero-locking of l-\ splittings for large mode linewidths. All 
of these problems have a common origin: a parameter is quite 
sensitive to the realization noise and the minimum of the likeli- 
hood function because of the realization noise can take a value 
which cannot be reached in the parameter range (here 0-90°). In 
our case to overcome this locking it is sufficient to use a param- 
eter j = sin 2 i instead of ;. If at the minimum, j is larger than 1 
it means that the fit using ; would have locked at 90°. Of course 
in this case it is not possible to recover a meaningful inclination. 
Obviously the locking at 0° is less likely to happen according to 
Fig. 6 since the black region is confined far from the limit: in 
other words, the likelihood surface along i — 0° is a ridge and 
not a valley as it is along ; = 90°. As observations become longer 
there are more and more points to describe the mode profile and 
the region is then more confined around the real minimum and 
further away from the 90° limit, the locking is then less likely 
to happen unless the underlying inclination is close to the limit 
itself. 

7. Conclusion 

As pointed out previously (Gizon & Solanki 2003; Ballot et al. 
2006) the extra parameter ; introduces difficulties in fitting stel- 
lar spectra, compared to the solar case. The main problem orig- 
inates in the blending of components within a multiplet which 
strongly correlates the inclination with the rotational splitting, 
making it difficult to disentangle them. For slow rotators, the 
uncertainties of i and v s are huge (§3.2), because of this. The 
projected splitting v* appears to be less sensitive to this, thus 
has more moderate errors. Beyond problems induced by correla- 
tions, other issues such as the locking of i at 90° often occur (see 
§6). The correct derivations of mode heights and linewidths are 
also dependent on a correct recovering of the inclination. 

These issues vanish for stars rotating sufficiently rapidly - 
or stars with p modes having sufficiently long lifetimes -, if we 
disregard other issues which appear such as problems in mode 
identification or blending between different modes. Fixing i be- 
fore fitting allows us to get around these difficulties. However, 
due to the numerous correlations between the different param- 
eters (see §4), the value of ; must be carefully chosen to avoid 
systematics. 

Nevertheless, even though the value of i modifies the error 
expected for the central frequency vo of a mode (§3.1), it does 
not introduce any biases as long as the multiplet is symmetric 
enough (§5). Hence, vo should be recovered in any case, as well 
as certain derived quantities such as the mode power. 
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